Pure Glass in Finite Dimensions 
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Pure glass is defined as a tiiermodynamic phase in which typical equilibrium particle configurations 
have macroscopic overlaps with one of some special irregular configurations. By employing 128- types 
of artificial molecules, a pure glass model is constructed in the cubic lattice. 
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Introduction Structural glass transitions have been 
extensively studied from several viewpoints. Among var- 
ious approaches, the random first-order transition sce- 
nario (RFOT) can be a candidate theory for the mean- 
field description of supercooled liquids and glasses [l|-y . 
However, the precise nature of glass in finite dimensions 
has not been clarified. For example, it is a long standing 
question whether a thermodynamic transition to a pure 
glass phase exists or not [5]. Why is such an apparently 
simple question difficult? In laboratory and numerical 
experiments, true equilibrium properties in prospective 
glassy materials are hardly observed, because the equili- 
bration time for a glass phase is considerably longer than 
accessible time scales. Theoretically, a renormalization 
group analysis may be a suitable approach to solving the 
above simple question [6|, l7| , but the agreement between 
the two studies has not been obtained. In this Letter, an 
attempt is made to answer the question of the existence 
of pure glass in finite dimensions by adapting a different 
approach. For this purpose, a slightly complicated model 
is constructed first, and then theoretical arguments con- 
cerning the model are presented. 

The first step in this approach is to assume statisti- 
cal mechanical conditions for constructing a pure glass. 
It is believed that particles in the glass phase are frozen 
into some irregular configurations. This hypothesis is ex- 
pressed by the following condition: typical equilibrium 
particle configurations have macroscopic overlaps with 
one of some special irregular configurations (MOSIC). 
The condition provides a precise definition of pure glass 
in this Letter. It should be noted that pure glass pre- 
sented thus far are limited to models in a random graph 
Sl- lipl o r with long-range interactions in the Kac limit 
111 . Il2| . One of the main results in this Letter is the 
presentation of a pure glass model with short-range in- 
teractions in finite dimensions. 

At this point, one may ask whether pure glass can be 
regarded as an idealization of glassy materials in nature. 
Furthermore, one may be interested in knowing whether 
the mode coupling theory [13| and the RFOT provide 
a useful description of the dynamics and thermodynam- 
ics of pure glass. An immediate answer to these ques- 
tions is not available. Nevertheless, since pure glass is 
different from gas, liquid, crystal, and quasi-crystal, the 
study on pure glass might shed light on the understand- 



ing of glassy materials and stimulate theoretical studies 
on glasses from a different perspective. 

Model If an infinite series of local minimum config- 
urations (LMCs) in a model are understood theoreti- 
cally, the statistical behavior of the model may be conjec- 
tured on the basis of an energy landscape of LMCs. The 
proposed model, where 128-types of artificial molecules 
are considered, is constructed along with this basic con- 
cept. Each molecule type is represented by an inte- 
ger a in {0, 1, . . . , 127}. Given a, the binary expan- 
sion cr = X]/i.=i cr*^^^2^~-^ uniquely defines a*^^^ G {0,1}, 
1 < k < 7. Further, a "sufficiently irregular" binary- 
array (cr*^^^((T))^^Q is picked and fixed. (Such an ar- 
ray was generated by selecting the value of cr'^^^(cr) as 
or 1 with probability 1/2 for each a.) As shown 
in Fig. [TJ the molecule type a corresponds to a unit 
cube whose k-th vertex is marked only when a^^^ = 1. 
Molecules are interpreted as three-dimensional general- 
ization of Wang tiles [14, 15], where a mark configuration 
such as {a'^^^ , a*^^^ , a*^^^ , a*^^^ ) on each plaquette represents 
the " color" of the plaquette. 
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FIG. 1. Unit cube with marked vertices. In the figure on the 



left-hand side, a 



(fc) 



1 only for /c = 1, 4, 5, 7. The figure on 



the right-hand side represents the top view of the cube, where 
one of the four colors is assigned for each pair of (a^ \a^ -*) 
aligned in the vertical direction. For example, red is assigned 



for 



.(2) 
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It is assumed that only one molecule can occupy each 
site in the cubic lattice A = {i e N^|z = (n,i2,^3),l < 
ik ^ L}' The total number of sites in the lattice 
is A^ = L^. A molecule configuration is denoted by 
o" = {o-i)ieA- For any nearest neighbor pair of sites (i, j) 



that satisfies the relation j — i = e/e, where e^ is the 
unit vector in the k direction of the lattice, the interac- 
tion energy ]4(<j, a') is defined for ((j^,crj) = (cr, a') as 
follows: Vk{(J,(J^) = 1 when mark configurations in the 
adjacent plaquettes of unit cubes at sites i and j = i-\-ek 
are different; otherwise, Vk{(J^(j') takes either or —1 ir- 
regularly. Such an irregular function was generated by 
selecting the value or —1 with probability 1/2 for each 
case. The irregularity of V is necessary to introduce a 
rugged energy landscape of LMCs. An example of the 
interaction energy is shown in Fig. [2l The Hamiltonian 
of the proposed model is given by 



H{a) = Y,Vk{cT^.cJ,), 



(1) 



{iJ) 



where (z, j) represents a nearest neighbor pair of sites. 
Although some elements of the interaction energy are 
determined by using probabilities, a quenched disorder 
does not exist in the lattice. Indeed, H{a) is trans- 
lational invariant if the periodic boundary condition is 
assumed. Throughout this Letter, P represents inverse 
temperature. 




FIG. 2. Example of the interaction energy. For two cubes 
cr = 29 (left) and a' = 106 (right), Vi{a,a') = 1, V2i(T,a') = 
1, V3{a,a') = 1, V2{a\a) = 1, and Vsia^a) = 1. Vi{a\a) is 
either or —1, which is selected with probability 1/2. 

Analysis A molecule configuration without any mis- 
matches between the mark configurations in all adja- 
cent plaquettes is called a perfect matching configura- 
tion (PMC). The advantage of the proposed model is 
that all PMCs can be constructed using a simple itera- 

(k) 

tion rule. The rule is as follows. First, the values of a^^ ^ 
for k = 1, 2, 3, 4 in plane is = 1; for /c = 1, 2, 5, 6 in plane 
22 = 1; and for /c = 1, 3, 5, 7 in plane ii = 1 are randomly 
chosen. Then, the molecule type di for z = (1, 1, 1) is de- 
termined. Accordingly, the value of cr| ^ for z = (1, 1, 1) 

(7) 

is given and this value is equal to the value of a^ for 
i = (2, 1, 1). Thus, the molecule type a^ for i = (2, 1, 1) 
can be determined. By repeating this procedure, all the 
values of ai for z G A can be determined depending on 
the initial choice of mark configurations in the planes. 
Since all PMCs are generated using the iteration rule, 
there are 2^^^^^^^^ PMCs. 



An important property of PMCs is that they exhibit 
neither long range positional order nor internal symme- 
try breaking. This property is evident from Fig. [3l where 
the cube configuration in plane is = L/2 is shown for a 
randomly chosen PMC. This configuration is certainly 
irregular. One can also examine this irregularity by mea- 
suring some statistical quantities. Theoretically, the ir- 
regularity may be understood from the deterministic it- 
eration rule used for constructing PMCs. Although its 
rigorous proof is not obtained yet, one can demonstrate 
that the rule maintains the random nature of configura- 
tions in the planes zi = 1, Z2 = 1, and is = 1 (See Ref. 
|16j for a related argument.). 
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FIG. 3. Perfect matching configuration (L = 32). Molecule 
configuration in plane is = L/2 is shown as the color repre- 
sentation of the top view of cubes. (See Fig. [1] for the color 
representation.) 

Next, the energy landscape of H{a) is discussed. It can 
be easily confirmed that the replacement of any cube in a 
PMC always increases the energy. That is, all PMCs are 
LMCs. Moreover, LMCs other than PMCs exist. Sup- 
pose that a cube at site i in a PMC is altered such that 
either a^ or a^ is changed. Then, three mismatches of 
ai are generated in the configuration. Two of them can 
be eliminated by employing the iteration rule beginning 
from the cube next to a^. This leads to an LMC with 
one mismatch which cannot be eliminated by any finite 
number of steps of cube replacement in the thermody- 
namic limit. The repetition of this procedure can yield 
LMCs with k isolated mismatches which cannot be elim- 
inated by any finite number of steps of cube replacement 
in the thermodynamic limit, where k is an arbitrary fi- 
nite integer. A collection of all such LMCs is denoted by 
{cr'^)aeA^ where A is the index set. \A\ is not estimated, 
but from the construction method of LMCs, it is found 
that 1^1 c^ exp(5oA^), where sq is a positive constant. 

Since the energy density Ua = H{cr^)/N for a G ^ 
is obtained as the space average of interaction energies 
for an irregular configuration cr'^^ the distribution of u^ 
has a sharp peak with width of 0{lJ_\fN). Thus, as 
studied in the random energy model [17|, it is expected 
for sufficiently low temperature that the statistical mea- 
sure condensates onto configurations near a cr^ satisfying 
H{a^) c^ u^N^ where u^ is the minimum energy density 
of LMCs in the thermodynamic limit. This condensation 
phenomenon leads to MOSIC. " 



Numerical experiments The condensation phe- 
nomenon is explored by numerical experiments. The 
statistical quantities for rather small-size systems in 
the equilibrium state under free boundary conditions 
are calculated by employing the exchange Monte Carlo 
method 0. 

First, thermodynamic quantities for the proposed 
model are investigated. As an example, in the left side 
of Fig. m the change in the energy density u = (u) 
with inverse temperature (3 is shown for L = 8, 9, 10, 11, 
where u = H{a)/N. An abrupt drop is observed around 
a temperature for each L. Since —du/dj3 is equal to the 
intensity of energy fluctuations Xu = N{{u — u)'^)^ Xu 
is displayed in the right side of Fig. HI The peak value 
of Xn, denoted by Xn^^^ increases as Xn^^ — ^^''^ (^^^ 
Fig. [5j), which is faster than xJT^^ — ^^ expected for the 
first-order transition. Furthermore, two graphs of Xu for 
L = 11 and L = 10 in Fig. |4] satisfy a relation 



Xu^L^^'f{{p-l3c{L))L'/n. 



(2) 



where a = 1.0 and u = 1/4.7 c^ 0.21. /3c(L) is esti- 
mated as the value of /3 where Xn becomes maximum. 
From these results, it is certain that a thermodynamic 
transition occurs at a non-zero temperature in the ther- 
modynamic limit, and it may be conjectured that the 
transition is the first order. It should be noted that the 
system sizes in the experiment are too small to identify 
the nature of the transition precisely. 




FIG. 4. (Color online) Thermodynamic quantities for several 
values of L. The left figure shows the energy density u as 
functions of /3 for L = 8 (asterisk), 9 (square), 10 (triangle) 
and 11 (circle). The intensity of energy fluctuations is dis- 
played in the right side. L = 9 (square), 10 (triangle) and 11 
(circle) . 

Next, in order to observe MO SIC in low tempera- 
tures, two independent identical systems are prepared, 
where two configurations are denoted by {a^ )ieA and 
{o'l )ieA^ respectively. The overlap between the two con- 
figurations is defined as 



ieA 



H<^1 



''\<rr') 



(3) 



Let P(g, /3) be the distribution function of q = q in the 
equilibrium state of the system with inverse temperature 
/3. In Fig. [6l the color representation of log P{q^/3) is 




FIG. 5. (Color online) The left graph indicates the maximum 
value of Xu for L. The guide line represents Xu = (L/3.4)^"^. 



The right figure shows XuL 
graphs with L = 10 and L 



~^'^ versus (/3- 
= 11. 



-/3c(L))L^-^fortwo 



shown, which clearly indicates the existence of two peaks 
at low temperatures. Since typical configurations in low- 
temperatures are irregular, it is concluded that MOSIC 
emerges. This is the most remarkable result in this Let- 
ter. 
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FIG. 6. A color representation of logP as a function of q for 
different values of /3. L = 10. 

A further evidence of MOSIC is proposed by study- 
ing the system under a special boundary configuration, 
where molecules in the boundary planes i^ = 1 and 
i/c = L (/c = 1, 2, 3) are fixed as those of a'' chosen from 
equilibrium configurations with large /3. Then, as shown 
in the left side of Fig. [71 the overlap with cr'' ^ which is 
characterized by 






(4) 



approaches one as /3 is increased, where A represents the 
inner region except for the boundaries in A. TV = (L — 2)^ 
is the number of sites in A. The result indicates that 
the boundary configuration selects one irregular con- 
figuration, just as the spin- up boundary configuration 
determines the ordered phase with positive magnetiza- 
tion. This is a direct evidence that there exists a low- 
temperature phase with MOSIC. 

Finally, properties near the transition point are briefly 
studied from a viewpoint of dynamics. The simplest 



quantity characterizing dynamical behaviors is the cor- 
relation function 



ieA 



(5) 



where cr^(O) is sampled in the equilibrium state and the 
time evolution a{t) obeys the Glauber dynamics with the 
heat bath method. It is observed that Cq{t) decreases 
exponentially in sufficiently high temperature, while it 
saturates to a finite value in the low temperature phase. 
Near the transition point, two types of trajectories, which 
correspond to de-correlation and freezing, coexist. The 
behavior may be consistent with the discontinuous tran- 
sition of q. Here, as shown in the right side of Fig. [71 
no intermediate plateau is observed, in contrast to the 
RFOT scenario. Although the decay of Cq{t) becomes 
slower as /3 approaches the transition point on the high 
temperature side, the divergent behavior is not concluded 
in the experiment of the small size system. 
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FIG. 7. (Color online) The left-side figure shows q* as a func- 
tion of /3 (L = 10). In the right side, Cq{t) are plotted with 
a log- log scale. /3 = 1.15 (squares), 1.155 (upward triangles), 
1.16 (downward triangles), 1.165 (diamonds), and 1.17 (cir- 
cles). Since trajectories for /3 = 1.17 are clearly classified into 
two groups, two graphs corresponding to the two groups are 
plotted. The dotted straight line represents a power- law form 
Cq(t) = 0.7t~°-'^ as a guide line. 

Concluding Remarks In this Letter, a pure glass 
model in finite dimensions has been proposed. In this 
model, typical equilibrium configurations in low temper- 
ature have macroscopic overlaps with one of some spe- 
cial irregular configurations. Although the behavior is 
similar to that observed in the models that exhibits one 
step replica symmetry breaking, the nature of the tran- 
sition seems to be different from the RFOT. In contrast 
to previous numerical studies on glass transition in finite 
dimensions [19-21], the present work has explored the na- 
ture of the low temperature phase more than properties 
near the transition point. The two peaks in P{q) in low 
temperatures may be the clearest demonstration among 
all existing studies of glassy systems in finite dimensions. 
Furthermore, the construction of an infinite series of local 
minimum configurations will provide a clue for theoret- 
ical analysis of statistical mechanics in sufficiently low 
temperatures. 

Before ending the Letter, I provide a few remarks. The 
arguments presented in this Letter are based on several 



conjectures whose validity is demonstrated with the aid 
of numerical experiments. A complete theory for describ- 
ing pure glass will be developed in future work. Then, 
it may be significant to construct a simpler model that 
exhibits a behavior similar to that of pure glass. Further- 
more, it is interesting to consider an effective field theory 
for describing pure glass, which might be related to the 
replica field theory [22]. 

The ultimate goal of this study is to construct a pure 
glass experimentally. Toward this goal, the next subject 
is to study a Hamiltonian system for which the rugged en- 
ergy landscape may be discussed theoretically. It will be 
amazing to observe the properties of pure glass via molec- 
ular dynamics simulations. In this study, complicated 
shapes of molecules will be employed so as to produce 
interaction potentials similar to 14 (a, a'). This challenge 
may be related to recent studies on complex structures 
emerging out of designed building blocks with anisotropy 
[23|, |2^. Here, let us recall a history of quasi-crystals. 
The mathematical study on aperiodic but regular tiling 
was conducted in the 1960s [14], whereas the experimen- 
tal construction of quasi-crystals was ffist reported in 
1984 [25]. As an extension of this Letter, a theoretical 
study on irregular but ordered tiling will be conducted. I 
hope that future researchers will attempt to consider the 
possibility of constructing pure glass experimentally. 
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